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Abstract 

The augmented space formalism coupled with the recursion method and a tight-binding 
linear Muffin-tin orbitals basis has been applied to study the effects of roughness on the 
properties of (001) surfaces of body-centered cubic Fe and face-centered cubic Co and Ni. 
The formalism is also proposed for the study of smooth surface. Comparisons have been 
made for three types of surfaces: a smooth surface, the surface with a rough top layer, 
and a more realistic model with several rough top layers converging into a crystalline 
bulk. Comparisons have been made between the magnetic moments, work function and 
electronic density of states in the three models described above. 

Keywords: Surface Magnetic Properties, Effect of surface roughness 


1. Introduction 

When a him is deposited layer by layer during an epitaxial growth, there is always a 
possibility of vacancies being created in the topmost layers. There are extensive studies 
on surface magnetism of ideal surfaces for over three decades. But a realistic surface 
is actually rough which forms during the process of epitaxial growth. Therefore, it is 
important to model a rough surface in any theoretical studies in realistic case. The 
combined effects of randomness in the incoming atomic hux and the smoothening due to 
surface diffusion in hnally creating rough surfaces and overlayers have been studied earlier 
[l|, 1^ . The top layers are expected to have maximum roughness as the atoms in these layers 
have looser bonding, lower coordination number and reduced symmetry as compared to 
bulk atoms. Rough surfaces can be modelled as a random binary alloys. We can envisage 
the missing atoms at a fraction of sites in the surface layers as random occupation of 
these sites by “empty spheres”. These spheres are hardly “empty”. They lack ion- 
cores, but carry non-negligible electronic charge. In the LMTO technique such empty 
spheres are introduced to take care of large interstitial space in non-compact structures. 
Consequently, when we solve the Kohn-Sham equations within the local spin-density 
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functional approximation (LSDA), the empty spheres are an integral part of the self- 
consistent treatment. We model our rough surface as a substitutional binary disordered 
alloy of the surface atoms carrying charge and ion-cores and empty spheres with charge 
alone. Due to these roughness, there is signihcant change in physical properties of a 
surface as compared to the bulk. 

We have taken two layers of empty spheres on top of the surface layer. This was to 
take care of charge leakage into the vacuum. Because of reduced co-ordination number 
at the surface, surface atoms are weakly bound compared to atoms in the bulk and so 
the interatomic separation is expected to be slightly different than in the bulk. Therefore 
lattice relaxation is necessary. In this communication we consider lattice relaxation of 
surface layers. A different relaxation procedure was adopted for Fe(OOl) where interlayer 
gap is relaxed using FPLAPW [^. 

To understand surface properties, it is very important to know the surface morphology. 
In transition metals, itinerant Fermi electrons play an important role in determining 
surface properties. In such systems, magnetic ordering arises due to the interaction among 
these Fermi electrons. Surface electronic and magnetic properties of such systems depend 
on the detailed electronic structure, the atomic arrangement and the composition. 

Our primary aim here to model realistic rough surface using Augmented space for¬ 
malism (ASF) but we have also shown that ASF can be extended to include ideal smooth 
surface. Therefore we have also carried out studies on ideal surface for comparison purpose 
with existing studies. 

2. Computational Methods 

Disorder has been tackled earlier with both the crude virtual crystal approximation 
and for a more sophisticated and self-consistent mean-field theory approach - the coherent 
potential approximation (CPA). The drawbacks of the CPA were immediately recognized 
and the search was on for going beyond mean-fields. The ASF jl, emerged as one of the 
very few methods that accurately described not only the averaged background disorder, 
but also a large deal of configuration fluctuations. The ASF has been described in detail 
in a number of papers (3,0, IJ and book j^. We refer the reader to them and mention here 
only those aspects of it which are relevant to our problem. In a nutshell, the ASR replaces 
the random parameters of the Hamiltonian by operators whose eigenvalues were the values 
taken by the parameters with probability densities which are the spectral densities of these 
operators. The augmented space theoremthen showed that the configuration average 
a special matrix element in configuration space. 



The right-hand side is the representation of an operator M in the configuration space 
$. The augmented Hamiltonian is 




RL 


RL 




RL,R'L' 


® {Srl,r'l'Trl,r'l' ® I) ® {^^Pr'l'® I— ^^)Pr'l'® Nr'^ (1) 


2 


Here P and T are projection and transfer operators. R are the lattice sites and 
L = {Im) are the orbital indices. For transition metal I < 2. and are the 

TB-LMTO potential parameters of the constituents atoms A and B of the alloy. S is the 
structure constant. 

The augmented space theorem gives the conhguration averaged Green function as 


< Grl,r'l'{z) >= {RL 0 0|(^/ - H)-^\R'L' 0 0) (2) 

The Green function is then obtained as a continued function expansion through the 
Recur sionMethod of Haydock et.al. and is generalized for operators on the real space 
augmented with the space of configuration fluctuations. It involves the familiar three 
term recursion : 


|1) = |R00) |2) = .H|l)-ai|l) 

|n + 1) = H\n) — an\Ti) — /S^_i|n — 1) 
an = {n\H\n)/{n\n) = {n\n) / {n — l\n — 1) (3) 


Leading to : 


^ Grr{z) ^ 




O'! 




z - a 2 - 




z-as- 




T{z) 


(4) 


a and (3 are potential parameters and T{z) is the terminator. 

The potential parameters are generated from TB-LMTO within local spin density 
approximation (LSDA) using Barth and Hedin exchange correlation potential. We have 
used seven shell augmented space calculation and nine steps of recursion. 


3. Results for Rough Surfaces 

We have chosen twelve atomic layers. The top most layer is relaxed by the 5% in¬ 
crement of lattice constant in the case of bcc Fe(OOl) whereas it is respectively 16% and 
9% in the case of fee Go(OOl) and fee Ni(OOl). We have considered two types of surfaces, 
one by roughening the top most layer and the other by roughening the top four layers 
with different degree of roughness. This hgure of percentage of increment are obtained 
by minimizing the total energy. 

Roughening is considered from 20% to 5% from top layer to the fourth lower layer 
with a difference of 5%. Figured] shows as we go down from the layer S to the bulk, width 
of density of states (DOS) increases. In the case of surfaces the width of the DOS of the 
top most surface layer (S) is narrower as compared to the bulk which is expected. The 
width increases substantially at S-1 layer and then reaches to the bulk values very slowly. 
Though width of DOS from S-2 layer onwards changes slowly, but there are significant 
changes in structure of DOS. This is due to variation in roughness. DOS bulk value is 
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Figure 1: Spin resolved density of states starting from the surface layer down to the bulk along the (001) 
direction. The roughness also decreases from 20% on the surface to none in the bulk. Fermi level is reset 
at zero. 4 




















































































Table 1: Work Functions. 


Properties 

Four layered rough surface 

Smooth surface 

Fe(OOl) 

Co(OOl) 

Ni(OOl) 

Fe(OOl) 

Co(OOl) 

Ni(OOl) 

Work Function (eV) 

4.15 

5.42 

4.79 

4.15 

5.33 

4.79 


Table 2: Layered based and bulk (B) orbital resolved magnetic moment in /is/atom. 


Layers 

Roughness 

(%) 

Fe(OOl) 

Co(OOl) 

Ni(OOl) 

s 

P 

d 

Total 

s 

P 

d 

Total 

s 

P 

d 

Total 

S 

20 

-0.06 

0.03 

2.82 

2.79 

-0.01 

0.03 

2.30 

2.32 

0.0 

0.0 

0.63 

0.63 

S-1 

15 

-0.05 

-0.03 

2.15 

2.07 

-0.07 

-0.01 

1.20 

1.12 

-0.02 

0.0 

0.43 

0.41 

S-2 

10 

-0.02 

-0.03 

1.89 

1.84 

-0.02 

-0.04 

1.82 

1.76 

-0.01 

-0.01 

0.64 

0.62 

S-3 

5 

-0.02 

-0.05 

2.13 

2.06 

-0.01 

-0.05 

1.97 

1.91 

-0.01 

-0.02 

0.74 

0.71 

S-4 

0 

-0.01 

0.0 

2.02 

2.01 

0.04 

-0.06 

1.72 

1.70 

0.01 

-0.02 

0.36 

0.35 

B 

0 

-0.02 

-0.06 

2.25 

2.17 

-0.02 

-0.06 

1.62 

1.54 

-0.01 

-0.02 

0.55 

0.52 


reached at the S-9*^ layer down the top most layer in the case of Fe(OOl) whereas it is 
reached at the S-8*^ layer in case of Co(OOl) and Ni(OOl). 

Figure [T] also shows interesting variation in spin resolved DOS. We note that down 
spin electrons contribute significantly to DOS at the Fermi level at the layers S & S-1 in 
the case of Fe(OOl). This trend is changed at the S-2 level where both up and down spins 
have significant contributions indicating that this layer is less magnetic. It is other way 
round at S-3 level where up has significant contribution and down is negligible. In the 
case of Co(OOl) contributions of spin-down to the DOS at the Fermi level is signihcant 
compared to spin-up at all the layers. This is also true in the case of Ni(OOl) except 
at S-1 layer where both have significant contributions. This means S-1 layer of Ni is 
almost non-magnetic at Fermi level. The effect of roughness is visible in this figure. For 
example in the case of Fe the width of DOS at S-4 layer where roughness is zero is more 
compared to other layer where roughness increases. Apart from the difference in width, 
the roughness smoothen the curve. 

Table [T] shows the work function. There is negligible effect of roughness on work 
function for the case of Fe and Ni but there is slight change in the case of Co. 

Table [2] shows that the topmost layer (S) of Fe(OOl) & Co(OOl) has the maximum 
magnetic moment. Whereas layer S-3 of Ni(OOl) has maximum magnetic moment. This 
is indeed true because splitting of spin up and down DOS is maximum at these respective 
layers as shown in the Figured] The bulk magnetic moment is attained at the S-9*^ layer 
in the case of bcc Fe(OOl) whereas it is attained at the S-8*^ layer in the case fee Co(OOl) 
and fee Ni(OOl). The average magnetic moment of the top four layers having different 
amount of roughness is more than its bulk value for all the three systems. This is expected 
as surface magnetic moment is enhanced compared to the bulk. The table [2] shows d- 
orbital alone contributes most to the the magnetic moment as expected. Similarly, Figure 
djshows d-orbital DOS is maximum whereas others are almost negligible. Among the three 
systems, the most significant change in d-band DOS of S layer is in Co(OOl). Therefore 
its magnetic moment in the S*^ layer is almost double that of S-1 layer. 

Majority of the work on surface properties is the study of rough topmost layer due to 
poisoning of substrate atoms 0,[I3, not due to voids. Therefore, we had roughened the 
top most layer with different degrees of roughness, (with 10% & 20% of empty spheres). 
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Figure 2: Orbital resolved total DOS for surface (S) and sub-surface (S-1) layers. Fermi level is reset at 
zero. 


Figure [3] shows that the peak near the Fermi level changes signihcantly when roughness is 
varied in all the three systems. Hence the electrons near the Fermi level are more affected 
by the change in roughness. It is found that the average magnetic moment decreases for 
all the systems when roughness is increased. When the roughness is more than 50% the 
surface magnetic moment increases and it approaches to the atomic magnetic moment 
for these systems. It is also observed that work function of the top layer changes slightly 
when roughness is varied. This value remains almost equivalent to smooth surface values. 
In this case also we observe the magnetic moment of the top most layer more than all 
other layers for Fe(OOl) and Co(OOl) where as layer S-3 has more magnetic moment than 
other layers for Ni(OOl). 

Comparison between the DOS of the two cases, i.e., a realistic case when first four 
layers are roughened (Figure 1) and the other when only the top layer is roughened with 
20% (Figure 3) indicates that there is minor difference in DOS in the case of Fe(OOl) and 
Ni(OOl). But there is significant difference in the case of Co(OOl). Similarly comparison 
of Tables |2] & |3] show magnetic moments get enhanced significantly for the realistic case 
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Figure 3: DOS of top layer (S) of Fe(001, Co(OOl) and Ni(OOl) in a model in which only this layer has 
been roughened only. Fermi level is reset at zero. 
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Table 3: Layer based and bulk orbital resolved magnetic moment when only top layer is roughened with 
20% in /xs/atom. 


Layers 

Fe(OOl) 

Co(OOl) 

Ni(OOl) 

s 

P 

d 

Total 

s 

P 

d 

Total 

s 

P 

d 

Total 

S 

0.01 

0.03 

2.81 

2.85 

-0.01 

0.03 

1.93 

1.95 

-0.01 

0.0 

0.65 

0.64 

S-1 

0.0 

0.0 

1.27 

1.27 

0.02 

-0.03 

0.92 

0.91 

0.02 

-0.01 

0.29 

0.30 

S-2 

-0.03 

-0.05 

1.82 

1.74 

-0.07 

-0.05 

1.75 

1.63 

-0.02 

-0.05 

0.47 

0.40 

S-3 

-0.03 

-0.04 

2.33 

2.26 

-0.02 

-0.05 

1.74 

1.67 

0.01 

-0.03 

0.76 

0.74 

S-4 

-0.03 

-0.05 

2.64 

2.56 

-0.01 

-0.06 

1.65 

1.58 

0.0 

-0.02 

0.55 

0.53 

B 

-0.02 

-0.06 

2.25 

2.17 

-0.02 

-0.06 

1.62 

1.54 

-0.01 

-0.02 

0.55 

0.52 


than for the other in the case of Co(OOl), where as there is only minor difference in the 
other two systems. 


4. Transition to smooth surfaces 


One would have thought, given the way we modelled the roughness of surfaces, that if 
we simply let the concentration of the empty spheres go to zero and we would recover the 
smooth surfaces. However, that is not the case. As the concentration of the empty sphere 
decreases, these ‘impurities’ become more and more isolated and form highly spikey im¬ 
purity states. The coherent potential approximation for example fails in this composition 
range and do not adequately reproduce the impurity structures in the density of states. 
Originally it was also thought that the ASR too misses out these structures. However, 
careful analysis of the “terminator” or the asymptotic behaviour of the continued fraction, 
indicated that this reproduces impurity peaks quite accurately. We have to incorporate 
not only the singularities at the band edges, but also those lying on the compact spectrum 
of H. Viswanath and Muller has proposed a terminator: 


T{z) = 


27r(E^)(P+2''+b/2 


B 


p + l 


A + q 


z-Eor{A-Ei){E2-z)y 


(5) 


The spectral bounds are at Ei, E 2 with square-root singularities, E:^ = E 1 E 2 and there 
is a cusp singularity aX Eq ii p = 1, q = 1 01 infra-red divergence ifp = —l/2,g = 0. Eq 


sits on the compact spectrum of H. Magnus [13[ has cited a closed form of the convergent 


continued fraction coefficients of the terminator : 


,2 _ 4n(n + q) _ 

(4n -I- 2g -I- p — 1) (4n -|- 2g -|- p -|- 1) 
,2 {2n + 2p + l){2n + 2q + p + 1) 

{An -I- 2^ -I- p -b 1) (4?7, + 2q + p + 3) 


( 6 ) 


The parameters of the terminator are estimated from the asymptotic part of the 
continued fraction coefficients calculated from our recursion. We shall use the Viswanath- 
Miiller termination appropriate for infra-red divergences and seamlessly enmeshed with 
the calculated coefficients as shown in Figure 01 


































Figure 4: (Top) Recursion coefficients from the recursive calculations (blue) and terminator (red) 
smoothly enmeshed. (Bottom) Density of states with a peak at the origin 

5. Discussions on almost Smooth Surfaces 
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Figure 5: Comparison of total density of states for bulk and surface (top most layer). Fermi level is reset 
at zero. 


Figure [5] shows energy bands of the surface states near Fermi level are narrower com¬ 
pared to bulk as expected. This is due to weakening of interaction by symmetry breaking 
and reduction in co-ordination number. Similar picture arise in the spin resolved DOS 
as shown in the Figure O Apart from narrowing of band there is change in number of 
spin-up surface states compared to the bulk states at the Fermi level. The amount of 
change is maximum in the case of Fe(OOl) and negligible in the case of Ni(OOl). But 
in all the three cases, there is signihcant change in the spin-down states. The splitting 
of spin-up and spin-down states near Fermi level is maximum for Fe(OOl) and least for 
Ni(OOl). This is expected because magnetic moment of Fe is maximum and it is least for 
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Ni. The magnetic moment is directly related to the amount of splitting in spin up and 
spin down sates. 
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Figure 6: Comparison of spin resolved density of states for bulk and surface (top most layer). Fermi level 
is reset at zero. 





Figure 7: Comparison of layer based total DOS for surface (S), sub-surfaces (S-1, S-2, S-3, S-4, S-5) and 
bulk (B). Fermi level is reset at zero. 


The layer based total DOS, shown in the Figure [7] indicates that DOS approaches 
the bulk value at S-5*^ layer in the case of Fe(OOl), whereas, it approaches to the bulk 
at S-T^ layer in the other two, Co(OOl) and Ni(OOl), cases. On comparison of figures [U 
and [7] we observe that the variation in structure in DOS between layers near the bulk is 
negligible in the case of smooth surface. The appearance of sharp peaks in DOS in figure 
[T] compared to figure [7] are due to disorderedness. For example, in the case of Co(OOl) 
20% produces very sharp spin-up peak compared to smooth surface. 

In transition metals, d-electrons contributes mainly to the formation of bands. There¬ 
fore, any surface effects must be reflected into d-band. This is evident from the Figur^ 
orbital resolved DOS. Apart from narrowing down of d-band, the band also gets splitted 
into few more prominent peaks at the surface .The narrowing of d-band is due to the 
dehybridization of s-, p-, and d-electrons fl3 |. 

Table lU shows magnetic moment approaches to the bulk value at S-5*^ layer in case of 
Fe(OOl) and at S-T^ layer in the other two cases as expected from DOS result. Magnetic 
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Figure 8: Orbital resolved DOS: surface : (a), (c) & (e) and bulk : (b), (d) & (f). Fermi level is reset at 
zero. 

moments for different layers exhibit Friedel oscillation. This is also shown in the figure [9l 
layer based percentage change in magnetic moment compared to the bulk value. Figure [9] 
further shows that the enhancement of surface magnetic moment is same for Fe(OOl) and 
Ni(OOl) but it is less for Co(OOl). Table 0] also shows that d-electrons contribute most 
to the magnetic moment compared to s and p electrons as expected form DOS. Figure [6] 
shows that at the Fermi level, difference between surface spin-up DOS and that of bulk 
is very small for Fe(OOl) and Co(OOl) and negligible for Ni(OOl). But there is significant 
difference between spin-down surface DOS and that of bulk DOS. This shows spin-down 
states are mainly responsible for the enhancement of magnetic moment at the surface. 
In all the three cases, the contribution of s- and p-orbital is comparable to TB-LMTO 
Green’s function method fl^ . 

For comparison we have summarized our results of magnetic moments along with 


Table 4: Orbital resolved magnetic moment for surface (S), sub-surfaces (S-1, S-2, S-3, S-4) and bulk (B) 
in /atom. 


Layers 

Fe(OOl) 

Co(OOl) 

Ni(OOl) 

s 

P 

d 

Total 

s 

P 

d 

Total 

s 

p 

d 

Total 

S 

0,01 

0.02 

2.80 

2.83 

-0.01 

0.02 

1.82 

1.83 

-0.01 

0.0 

0.69 

0.68 

S-1 

-0,02 

-0.03 

2.13 

2.08 

-0.01 

-0.03 

1.49 

1.45 

-0.01 

-0.01 

0.60 

0.58 

S-2 

-0,03 

-0.06 

2.06 

1.97 

-0.01 

-0.04 

1.79 

1.74 

0.0 

-0.02 

0.62 

0.60 

S-3 

-0,02 

-0.04 

2.70 

2.64 

-0.01 

-0.06 

1.66 

1.59 

0.0 

-0.02 

0.56 

0.54 

S-4 

-0,01 

-0.06 

2.43 

2.36 

-0.02 

-0.06 

1.62 

1.54 

-0.01 

-0.02 

0.55 

0.52 

S-5/B 

-0,02 

-0.06 

2.25 

2.17 
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Figure 9: Percentage variation of magnetic moment with respect to the bulk value for different layers. 


other available calculated and experimental values in table O Table E] shows that the 
enhancement of surface magnetic moment compared to the bulk is 30% for Fe(OOl), 19% 
for Co(OOl) and 32% for Ni(OOl). Our result of bulk magnetic moment for Fe is within 
5% difference from experimental value B Bl- Since enhancement of surface magnetic 
moment is less for Co(OOl) than the other two, therefore it is more stable to the change 
in the environment. 

Experimental values of local magnetic moment, like surface layer and layers below are 
not available for comparison. It is only the bulk magnetic moment for which experimen¬ 
tal result is available. Work function is a surface property and is related to the surface 
DOS. It is an experimentally measurable quantity and therefore theoretical calculation of 
it is necessary for comparison with experimental result to test the accuracy of theoretical 
study. Table [H] shows our calculated work functions for all the systems under consideration 
agree quite well with experimental and other theoretical calculations. 


6. Conclusion 

The augmented space formalism coupled with the recursion method and TB-LMTO 
has been successfully applied to study the effects of roughness on surface properties of 
bcc Fe(OOl), fee Co(OOl) and fee Ni(OOl). Two types of roughening are considered here, 
one with 10% & 20% roughening only at the top layer and the other a more realistic 
surface, i.e., by roughening the hrst four layers with 20%, 15%, 10%, & 5% from top 
layer respectively. Twelve atomic layers are considered to fully achieve bulk properties. 
Two layers of empty spheres are considered above the surface to include the vacuum. It 
is found that the trend in the variation of the width of DOS among the layers changes 
when a realistic surface is considered in comparison to a smooth surface. The comparison 
of surface properties for two different types of roughed surfaces shows slight difference 
in surface magnetic moment of Fe(OOl), whereas signiheant difference is obtained for 
Co(OOl) and no change is in the case of Ni(OOl) between these two types of surfaces. The 
magnetic moment at the top layer is maximum in the cases of Fe(OOl) and Co(001) but 
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Table 5: Comparision of magnetic moment in for surface (S), sub-surfaces (S-1, S-2, S-3) and 

central layer or bulk (C/B). Number in the square brackets represents the reference numbers. 


Methods 

Fe(OOl) 

Go(OOl) 

Ni(OOl) 

S 

S-1 

S-2 

G/B 

S 

S-1 

S-2 

S-3 

S-4/B 

S 

S-1 

S-2 

S-3 

S-4/B 

FP- 

LAPW 

2.98 

[18,i9, 20] 

2.35 

\m 

2.39 

\m 

2.20 

m 

1.86 

mi 

1.64 

[2^ 

1.65 

[21J 

1.64 

[2^ 

1.65 

mi 

0.68 

\m 



0.56 

m 


2.80 

[3] 

2.38 

[20] 

2.43 

[20] 

2.15 

m 






0.73 

[20] 

0.68 

[20] 

0.66 

[20] 

0.63 

[20] 





2.30 

[20] 






0.68 

[22] 

0.60 

[22] 

0.59 

[22] 

0.56 

[22] 


LMTO 

2.87 

]23,2£ 

2.34 
[23, 2^ 

2.33 
[23, 2^ 

2.18 
[23, 2^ 






0.59 
[23, 2^ 

0.58 
[23, 2^ 

0.57 
[23, 2^ 

0.55 
[23, 2^ 


TB- 

LMTO 

2.97 

[la 

2.30 

\m 

2.37 

\m 

2.24 

m 

1.84 

[la 

1.63 

\m 

1.66 

[1^ 

1.64 

[la 


0.69 

[la 

0.64 

\m 

0.66 

\m 

0.64 

[1^ 


2.97 

2.30 

2.37 

2.25 

[2^ 

1.84 

[2^ 

1.63 

1.66 

[2^ 

1.65 

[2^ 

1.66 

[2^ 

0.69 

0.64 

0.66 

0.64 

0.65 

[2^ 

2.86 

[9] 

2.16 

[9] 

2.38 

[9] 

2.17 

[9] 

1.76 

m 

1.46 

m 

1.58 

[26] 

1.56 

[26] 

1.58 

[2^ 

0.65 

[26] 

0.53 

[26] 

0.61 

[26] 

0.60 

[26] 

0.59 

[26] 

2.98 

[27] 

2.17 

[27] 

2.40 

[27] 

2.26 

[27] 











2.95 

m 

2.20 

m 

2.39 

\m 

2.28 

\m 











2.99 

[27] 

2.21 

[27] 

2.38 

[27] 

2.26 

[27] 











LCAO 

3.01 

[2^ 

1.69 

m 

2.13 

m 

1.84 

m 






0.44 

[29] 

0.58 

[29] 

0.62 

[29] 

0.56 

[29] 

0.54 

[29] 

ASR 

2.99 

[9] 

2.17 

[9] 

2.38 

[9] 

2.27 

[9] 











EXPT 




2.21 

m 





1.71 

[17] 





0.616 

m 




2.22 

[17] 











OUR 

WORK 

2.83 

2.08 

1.97 

2.17 

1.83 

1.45 

1.74 

1.59 

1.54 

0.68 

0.58 

0.60 

0.54 

0.52 



Note that Refs [2 


2bll use 

dilatation; Refs 121 


Green’s function techniques; Ref 
are Supercell calculation; Refs 127 


is with TB-LMTO including surface 
l26l | use Real space recursion. 


Table 6: Work functions of Fe(OOl), CofOOl) and Ni(OOl) in eV. Number in the square brackets represents 
the reference numbers. Note that Ref [l^ is with Green’s function; Refs 00,133 use surface embedded 


Green’s function. 


Methods 

Fe (001) 

Co(OOl) 

Ni(OOl) 

FPLAPW 

4.29 [1^ 

5.17 [^ 

5.37 [22], 
5.5 M, 
5.31 [30] 

LAPW 



5.71 [31] 

LMTO 

4.30 [2^ 


5.02 [2^ 

TB-LMTO 

4.5 [1^ 

5.52 [1^ 

5.75 [1^ 

Experimental 

4.67 [32], 
4.4 [3^ 

5.0 [32] 

5.22 [32] 

Our work 

4.15 

5.33 

4.79 
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it is at the layer in the case of Ni(OOl) in both types of roughening. Layered based 
magnetic moment is also found different in both types of surfaces. Work functions of the 
respective systems found to be almost same for both types of rough surfaces as well as 
for the smooth surface. In all the cases, d-band has signihcant contribution towards the 
magnetic moment and there is signihcant surface effects on the d-band DOS in all the 
systems. We have shown that Augmented space formalism which works well for rough 
surfaces can also be applied to the nearly smooth surface. Therefore, we have carried 
out detailed study of the same systems for smooth surfaces also. Our result agree quite 
well with other available results. Layer wise magnetic moments show Friedel oscillations. 
The work function is found to be 4.08 eV, 5.30 eV and 4.76 eV for Fe(OOl), Co(OOl) 
and Ni(OOl) respectively. These values are closed to the available experimental and other 
calculated results. 
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